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Abstract 



We obtain the time-dependent correlation function describing the evolution of a single spin 
excitation state in a linear spin chain with isotropic nearest-neighbour XY coupling, where 
the Hamiltonian is related to the Jacobi matrix of a set of orthogonal polynomials. For the 
Krawtchouk polynomial case an arbitrary element of the correlation function is expressed in 
, a simple closed form. Its asymptotic limit corresponds to the Jacobi matrix of the Charlier 

' polynomial, and may be understood as a unitary evolution resulting from a Heisenberg group 

^ , , element. Correlation functions for Hamiltonians corresponding to Jacobi matrices for the Hahn, 

dual Hahn and Racah polynomials are also studied. For the Hahn polynomials we obtain the 
general correlation function, some of its special cases, and the limit related to the Meixner 
polynomials, where the sm(1, 1) algebra describes the underlying symmetry. For the cases of 
^ , dual Hahn and Racah polynomials the general expressions of the correlation functions contain 

summations which are not of hypergeometric type. Simplifications, however, occur in special 
cases. 



m 
p 

^ '. 1 Introduction 

o\ 

■ Transfer of a known or an unknown quantum state from one site to another is a key requirement in 
linking neighbouring small quantum processors for facilitating large scale quantum computation. 

k> I S. Bose [1,2] introduced linear spin chains as a channel for such short distance quantum commu- 

j_j ■ nication. Such a connector has an inherent advantage as it renders the quantum processor and 



the communicating channel to be made of the same physical system. This eliminates the need of 
developing interfaces. Transmission of data in such linear quantum registers has been the subject 
of many investigations [1,3-6]. Using the classical concept of group velocity of a wave packet it 
has been observed [7] that one-dimensional spin rings allow high-fidelity transmission of quantum 
states if both communicating parties have access to a finitely limited number of qubits in the ring. 
Quantum communication with closed spin ring with twisted boundary conditions has been dis- 
cussed in Ref. [2]. An excellent review [8] and references therein describe the current developments 
in this field. 

The transmission of quantum states can in principle be performed by a chain of qubits coupled 
via the Heisenberg or the XY interactions [9-12]. Interesting situations arise if one assumes to 
have individual control of the nearest-neighbour couplings in the spin chain. The idea of pre- 
engineered interqubit couplings has been discussed considerably [13,14]. One of the advantages of 
well-chosen controlled couplings is that one can obtain mirror inversion of a quantum state with 
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respect to the center of the chain, and that perfect transfer of quantum states is possible [3,6, 15] 
at certain specified times over arbitrary length of the spin chain. Propagation of entangled states 
in anisotropic XY spin chains has been discussed in [16]. Using a system based on a dispersive 
qubit-boson interaction to mimic XY coupling that relaxes the nearest-neighbour restriction the 
transfer fidelity of the chain has been found [5] to achieve a nearly optimal value. 

One of the main results of [15] is the introduction of two (analytic) mirror-periodic Hamiltonians 
that allow for perfect transfer. To achieve perfect mirror inversion of an arbitrary many excitation 
state it is sufficient to consider transformations of all single excitation states to their mirror images 
as the time-dependent transition amplitudes for multiple excitation states may be constructed [15, 
17] via the Slater determinant of their single excitation counterparts. The eigenstates of the single 
excitation sectors of these Hamiltonians are related [15] to discrete orthogonal polynomials, namely 
Krawtchouk polynomials and dual Hahn polynomials. 

In the present paper we investigate such systems from a general point of view. We shall assume 
that the chain of qubits is described by a Hamiltonian of XY type, in such a way that the interaction 
strengths (and qubit energies) are related to the Jacobi matrix of a system of discrete orthogonal 
polynomials. This will allow us to give a general expression for the transition amplitude of a single 
excitation from the sth site (sending site) to the rth site (receiving site) in the chain of spins. For 
some types of discrete orthogonal polynomials, the expression of these transition amplitudes can 
be simplified, leading to some interesting properties. 

Let us consider a by now classical system of + 1 interacting qubits (spin 1/2 particles) in a 
quantum register, with an isotropic Hamiltonian of XY type: 

JV-l N 
fc=0 k=0 

where Jk is the coupling strength between the qubits located at sites k and k + 1, and hk is the 
"Zeeman" energy of a qubit at site k. So the subindex k {k = 0, 1, 2, ... , A^) labels the position 
of the qubit in the chain, and the superindex refers to the Pauli matrices a^, and . The 
Hamiltonian ([1]) preserves the total spin: [H , X^^g ~ therefore, we may analyze 

various spin excitation sectors separately. 

To describe the Hilbert space associated with the Hamiltonian, one adopts a standard fermion- 
ization technique [17]. The Jordan- Wigner transformation [18] maps the Pauli matrices to spinless 
lattice fermions {a^, /c = 0, 1, . . . , N} obeying the anticommutation rules 

{4,0£} = 4,£, {afc,a£} = {4,4} = VA:,£g {0,1,--- ,iV}. (2) 
We may now recast the Hamiltonian ([1]) as 

N-l N 

H Jk{a\ak+i + al+i^fc) + E ^ka\^k (3) 

that describes a set of + 1 fermions on a chain with nearest-neighbour interaction (hopping 
between adjacent sites of the chain), and subject to a non- uniform background magnetic field 
denoted by (k = 0, 1, . . . , A'"). We shall assume that the system is initially in its completely 
polarized ground state |0) = |00---0) = |0) (8) |0) (55 • • • |0), where |0) denotes the spin down 
state. Let \k) = \00 ■ ■ ■ 010 • • • 0) = a|.|0) (A; = 0, 1, ... , A^) denote a state in which there is a single 
fermion at the site k and all other sites are empty, i.e. \k) describes the state in which the spin at 
the site k has been flipped to |1). Clearly, the set of states \k) (fc = 0, 1, . . . , A^) forms a basis for 
the single-fermion states of the system, and we may represent them by the standard unit vectors 
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in column matrix form: 







\k) 



0,1,... ,A^). 



(4) 



In this single- fermion basis, the Hamiltonian H takes the matrix form 





( ho 
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hi 


Ji 







M = 





Ji 


h2 














Jn-i 






V 







Jn-1 h^ 


/ 



(5) 



The dynamics (time evolution) of the system is completely determined by the eigenvalues ej and 
eigenvectors cpj of this matrix. It is then, as noted before, a standard technique [15, 17] to de- 
scribe the n-fermion eigenstates oi H {n < N) using the single-fermion eigenstates and Slater 
determinants. For this reason we concentrate here on the single-fermion eigenstates. 

The matrix M in ([5]) is real and symmetric, so the spectral theorem [19] implies that it can be 
written as 

M = UDU^ (6) 



where D is a diagonal matrix and U an orthogonal matrix: 

D = diag(eo,ei, . . . ,eiv) 



U^U = I. 



(7) 
(8) 



The entries of D are the single-fermion energy eigenvalues, and the columns of the matrix U are 
the (orthonormal) eigenvectors of M, i.e. the single-fermion eigenstates: 



Ui, 
V Un, J 



N 



N 



Ukj \k) = Ukj 4|0) (j = 0, 1, . . . , N), 



(9) 



k=0 



k=0 



with H(f>j = M(j)j = € 



v^j. From the orthogonality of U, the inverse relation follows: 

N 

\k) = J2^kj(t>j- 

j=0 



(10) 



We now turn to the dynamics of the system under consideration, described by the unitary time 
evolution operator 

L{{t) = exp{-itH). (11) 

Assume that the "state sender" is located at site s of the spin chain, and the "state receiver" at site 
r (s and r are site labels, belonging to {0, 1, . . . , N}). At time t = the sender turns the system 
into the single spin state \s). After a certain time t, the system evolves to the state U{t)\s) that 
may be expressed as a linear superposition of all the single spin states. So the transition amplitude 
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of an excitation from site s to site r of the spin chain is given by the time-dependent correlation 
function 

frAt) = {rmt)\s). (12) 

This is the central quantity of this paper (sometimes referred to as the "correlation function"), and 
the square of its modulus gives the transition probability from the sth to the rth spin excitation 
state. Note that it can be expressed by means of the orthogonal matrix U appearing in dH). Indeed, 
using the definition (jl2p . the expansion (jlOp . and the orthogonality of the states one finds: 

N N 
fr,s{t) = {y~^ Urk(t)k\ ex.p{-itH) Usj4>j) 

k=0 j=0 
N N 
= Urk(pk\ ^ UsjC 

k=0 j=0 

N 

= UrjUsje ^ . 
j=0 

In other words, using the abbreviation z = e~**, one has 



-itei 



(13) 



frAt) 



N 



(14) 



The purpose of this paper is to show that various interesting closed form expressions can be given for 
this crucial quantity /r,s(t), in the case that the fixed values characterizing the system (the values 
Jk and hk) are related to the Jacobi matrix of a set of discrete orthogonal polynomials. We shall 
illustrate this first by means of an example, where the polynomials are Krawtchouk polynomials. 
From this example, the general technique will be clear. Then we continue analyzing some systems 
related to other classes of orthogonal polynomials. 

From the outset it is assumed that such chains of qubits can be pre-engineered for any given 
set of values Jk and /i^, i.e. the strength of the couplings can be engineered and the external static 
potential at site k is controlled. Such a physical control has been the subject of many papers (see [3] 
and references therein). 

We end this section by two remarks. First of all, the strengths Jk in ([3]), and thus the off- 
diagonal elements of M in ([5]), are positive. From the mathematical point of view, however, the 
problem can just as well be solved for a matrix with negative off-diagonal elements. Indeed, if 



M' 



ho 


-Jo 








-Jo 


hi 


- Jl • • • 








-Jl 


h2 





V 







-Jn- 



\ 



-Jn-1 

hjy ) 



(15) 



then M'-^ = (-l)^+'=Mjfc. This implies that M' has the same eigenvalues as M. Moreover, 
the orthogonal matrix U' with U'-^ = {—^})^^^\Jjk diagonalizes M' in the same way as in (l6|), so 
the eigenvectors of M' (i.e. the columns of C/') are the same as those of M up to sign changes in 
the components. Then the transition amplitude follows from (jl4p and is essentially the same as 
that corresponding to M: /^^(t) = (— l)'"~'~^/r,s(i). Secondly, also a matrix that differs from the 
original one by a constant factor and a multiple of the identity matrix leads essentially to the same 
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computation. Indeed, let M' = AM + fil, where A and fi are constants and / is the identity matrix. 
Then the same matrix U from ([6]) diagonahzes M' , the only difference being the eigenvalues given 
now by Ae^ + fi. Also, it follows immediately from that frg{t) = e~**'^/r,s(At). 

2 The Jacobi matrix of Krawtchouk polynomials 
2.1 Computation of the general correlation function 

Let us start this section by introducing some standard notation and known facts of Krawtchouk 
polynomials. We follow the notation of [20]; other standard works on (discrete) orthogonal poly- 
nomials are [21, 22]. The Krawtchouk polynomial of degree n (n = 0, 1, . . . , N) in the variable x, 
with parameter < p < 1 is given by 



-X, —n 1 
-N 'p 



(16) 



The function 2-^1 is the classical hyper geometric series [23,24], and in this case it is a terminating 
series because of the appearance of the negative integer —n as a numerator parameter. Krawtchouk 
polynomials satisfy a (discrete) orthogonality relation [20]: 



N 



'^w{x)Knix)Kmix) = dnS 

n 



(17) 



x=0 



where w{x) is the weight function in x and d„ is a function depending on n: 

N-x n . An. J \ (\-V 



W[X} 



(1 - p) 



(x = 0,l,...,iV); 



dn 



They also satisfy the following three-term recurrence relation: 

- xKnix) = n(l - p) Kn-i{x) - [p{N - n) + n(l - p)] Kn{x) + p{N - n) K„+i(x). 
It is often convenient to introduce orthonormal Krawtchouk functions by 

^/w{x)Kn{x) 



Kn{x) 



(18) 



(19) 



(20) 



Then, rewriting the orthogonality relation and the recurrence relation in terms of the functions 
Kn{x), it is easy to obtain the following: 

Lemma 1 (see [25]) Let Mx be the tridiagonal (N -|- 1) x (N + l)-'matrix (Jacobi matrix) 

\ 





( ho 


-Jo 







-Jo 


hi 


-Jl 


Mk = 





-Jl 


h2 



where 





—Jn-1 
-Jn-1 hN J 



Jn = Vpi'^-p)Vin + l){N -n), hn = Np + {l- 2p)n, 



(21) 



(22) 
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and let U be the {N + 1) x (N + l)-matrix with elements Ujk = Kj{k). Then 

UU^ = U'^U = I and Mk = UDU^ (23) 

where 

I) = diag(0,l,2,...,A^). (24) 

In other words, the eigenvectors of the Hamiltonian (in the single-fermion case) corresponding 
to the quantities (|22p have components equal to normahzed Krawtchouk polynomials, and the 
corresponding energy eigenvalues are ej = j (j = 0, 1, . . . , N). Note that Krawtchouk polynomials 
have been used before as a basis for quantum chains [26]. Here, only their evaluations at integer 
values of the support are used as matrix elements of U in the diagonalization process. 
Let us now consider the transition amplitude or correlation function: 

N N 

k=0 k=0 
1 ^ 

So we need to compute the quantity in (j25|) . First of all, note that in general fr,s{^) is a periodic 
function of t since z = e~**. In particular, it follows from (j25p and the orthogonality relation (|17p 
that fr,s{t) = Srs for t = and for any multiple of 27r. So after a time span of 2tt, the system is 
back in its original state where only the spin at the sending site s is flipped. 

The purpose is now to compute (I25p explicitly. We shall do this in two ways: a classical method 
and a group theoretical method. The classical method is short and straightforward. Rewriting 
the polynomials in ()25p as 2-P'i-series, this sum reduces to a classical summation formula given for 
example in [27, p. 84, (8)]. This leads immediately to the following closed form expression: 

(26) 

Before discussing some special and interesting cases of this formula, we shall also deduce this 
correlation function in a different way. 



2.2 Group-theoretical computation 



The group-theoretical way to obtain (j26l) is somewhat longer, but it does not use any reference to 
orthogonal polynomials or summation formulas of hypergeometric type. So it sheds another light 
on why the final formula (j26p is so simple. 

Consider the Lie algebra su{2) of quantum angular momentum theory [28], with basis Lq, L± 
and commutation relations 



[Lo,L±]=±L±, [L+,L_]=2Lo. 
For any positive integer, the {N + l)-dimensional irreducible representation is given by 

\ 



(27) 



/ ^ 

/ 2 



TV _ 1 
2 ^ 



N 



+ 1 



N 
' 2 



/ Vl • 

^/2-iN- 1) 







(28) 
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and L_ = (L_|_) . Note that, in this representation, the matrix Mx from (|2ip is written as 

N 



Mk = -I + {2p - l)Lo - - P) {L+ + L-) , (29) 

where I is the identity matrix. Consequently, the computation of 

/,,,(*) = (r I exp(-itMi^) I s) (30) 

simply leads to the computation of a matrix element of an SU (2) group element. So apart from a 
factor e~**^/^, we need the computation of the following left hand side, which we write as 

exp(-if(2p- l)Lo + itVP(l^^(-^^+ + ^-)) = e«-^-e''^Oe^^+. (31) 

The equation (j3ip is an example of the standard Baker-Campbell-Hausdorff (BCH) type of de- 
composition of an SU{2) group element. For a general discussion on BCH decomposition of group 
elements see [29]. As we intend to use such decompositions in other contexts, we provide here 
the usual procedure of derivation of these rules. The constants ^, r/, C, in ([3T|) should be repre- 
sentation independent (all elements are group elements), so to determine these we can perform 
the computation using any faithful representation, and, in particular, the standard 2-dimensional 
representation. Expanding the exponential in the left hand side of (j3ip . using the 2x2 matrices 



one obtains after some calculations: 

cos| -i(2p- l)sin(l) 2i^p{l-p) sin(f) \ 

2i^p{l-p) sin(f) cos(f) +i(2p- l)sin(|) y ■ ^ ' 

On the other hand, the right hand side of (j3ip yields, in the same 2-dimensional representation (I32p : 
Identification yields: 



V2 (^\ 70 ^\ ■ f*\ c A 2i^yp(l-p) sin(|) 

3 " = cos - — z(2p — 1) sm - , t=C = TTs : ^ — , , , . (35) 

^2^ ^ ^2^' ^ ^ cos(|) -i(2p- l)sin(l) ^ ^ 



Now we compute the matrix element of the right hand side of (13ip for an arbitrary representation 
of dimension {N + 1). Using the common matrix elements (|28p . one finds 



Lo\s) = {--s)\s), ^|.) = . J \s-k), (36) 



2 ' ' kl ' ' \l \kj \ k 

and similarly 



(37) 
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Now (using ^ = C) 



^ Cj+k 

(r|e«^-e''^Oe^^+|s) = V ^-{r\{L_Y e''^^ {L+f\s) 
^-^ j\k\ 



j,k=0 



rlsl ^s-rn(Il-r)^ {N-r + j)\ 



{N - r)\{N - s)r y 3\{r - 3)l{j + s - r^}^ ' ' ^ ' 

The last sum is of hypergeometric type and is proportional to 

( -N'^^)='^' ( -iv'up(l-p)sin^(i)) ' ^^^^ 

using ([35]) . But z = e~**, so (1 — z)^ /z = — 4sin^(|), and taking all factors together one recov- 
ers p6p . Note that this computation is essentially equivalent to that relating matrix elements of 
finite rotations (d-functions) to Jacobi polynomials [28, Chapter 4]. 



2.3 Discussion and special cases 

Let us return to the closed form expression (j26]) for the transition amplitude fr,s{i) at time t. As 
already noticed, since fr,s(t) is a function of z = e~**, the system is periodic with period 27r. At the 
initial moment t = the system is in the state with all spins down except at site s where the spin 
is up. At any time which is a multiple of 2tt, the system is back in this initial state: fr^sC^kir) = 5rs 
for k £ Z4_. At other times, the system is in general in a mixed state. Due to the orthogonality of 
the basis states, one has 

AT 

r=0 

for any s and any time t. In fact, more generally, the matrix of correlation functions is unitary, so 
for any time t and any indices r and s: 

AT N 
k=0 k=0 

Let us consider the case when the sender is located at site 0, i.e. s = 0. Then (p6]) yields 

fr,o{t) = ^/(^(vW^r(i - zni-p+pzf-^. (40) 

So far, p (0 < p < 1) is still a free parameter. A special case occurs when p = 1/2: 

frAt) = (^) (1 - + 'f"'- (P = 1/2) (41) 
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Using z = e **, this gives 



N-r 

■ iP = 1/2) (42) 
In other words, 

/r-,o(vr) = <5,,iv. (p=l/2) (43) 

This is the situation of "perfect state transfer" described aheady in [15]: at time t = tt the system 
is in the state with all spins down except at site the spin is up. So for this time there is perfect 
state transfer from site to site N. 

Let us mention here that the condition for perfect state transfer can be deduced from the 
corresponding Jacobi matrix itself [6]. In order to allow perfect state transfer, the matrix ([5]) 
should be mirror-periodic, i.e. /i„ = hjy^n and J„ = Jn-i-u for all n. Clearly, for ([2T]) this is the 
case when p = 1/2, see (j22]) . 

More generally, let us specialize the expression (I26p for time t = tt: 

^■-<" = \/(^)(")(vW^r'^'-(i - ^pf-'-' {":-'■• 4^) . (") 

This expression shows once again that taking the free parameter p = 1/2 yields a special case: 

/r,s(7r) = 6r+s,N- {P = 1/2) 

So for p = 1/2 there is also perfect state transfer between the sites s and N — s. 




2.4 A limiting case 

A classical limit of Krawtchouk polynomials are Charlier polynomials. Putting the parameter 
p = a/N, and letting N go to +00 yields Charlier polynomials C„(x;a) [20]: 

lim Kn{x; -,N) = C„(x; a) = 2F0 { ' ; , (45) 



^^^^^ --/tV 'Ar" / --/tv— 7— / i^ui •> 

i\ \ — a 

satisfying the orthogonality relations 



00 ^ I 



y^—e'^Crn{x;a)Cn{x;a) = —5mn, (46) 

x=0 

where a > is a positive parameter. The recurrence relation reads 

xCn{x; a) = -aCn+i{x; a) + (n + a)Cn{x; a) - nC„-i(x; a). (47) 

To see which spin chain corresponds to this limit, we take the appropriate limit in ()22p . and 
find for 

H = ^(a + k)alak - ^ a{k + l)(4afc+i + al.+i«fc)- (48) 

k=0 k=0 

The above infinite chain of fermions with nearest-neighbour interaction immediately provides a 
unitary representation of the (central extension of the) Heisenberg algebra f)4(a, a^, DT, J), where 
the algebraic generators may be constructed as 

00 00 00 00 

g = Vfc + 1 alak+1, = Vfc + 1 Qfc^iQfc, DT = y^A;a^afc, J = ^a^afc. (49) 

k=0 k=0 k=0 k=0 
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Employing the defining anticommutation relations of the fermionic variables ^ we observe that 
the above operators satisfy the commutation relations of the Heisenberg algebra Il4(a, o^,DT, J): 

[o, a^"] = J, [% a] = -0, [% = [X, 3] = where X G {DT, o, a^}, (50) 

where the unitary representation in terms of the single spin excitation states (]10p reads 

\k) = VkTT\k + l), a\k) = ^\k-l), ^\k)=k\k). (51) 

Using the generators introduced in (I49p we may express the Hamiltonian (j48p as 

H = ^ + a3 - ^/aia + a^). (52) 

Following [29] the time evolution operator introduced in (jlip may now be expressed in the normal 
ordered BCH-factorized form: 

h({t) = exp{a{z — 1)) exp(^/a(l — z)a^) exp(-y/a(l — z)a). (53) 



It is worth mentioning that the correlation function (jl2p obtained via the operator factoriza- 
tion ([55]) . and the appropriate limiting value easily computed from ([25]) precisely agree. We quote 
the final answer: 



As(i) = \/^(l-^r^e-"+-2Fo( " (54) 



r+s 



r!s! \ — ' a(l — z)^ 

Intuitively the above result is easily understood from the well-known [29] result that the large spin 
(y) contraction limit of the su{2) algebra ()27p is given by the Heisenberg algebra. From ()54p it 
follows that the state with spin up at position s = emanates over the infinite chain with the 
correlation function given below 

/,0(t) = J^{1 - g-it)rg-a+aexp{-if)_ (55) 

Starting at time t = 0, some reflection takes place at the "infinite end" of the chain, and by time 
t = 2Tr the system is back in its original state. In the middle of this, at time t = n, the state is 
"spread" over the infinite chain according to the amplitude 



-2« ./(4«)^ 



/r,o(7r) = e 

V r! 

Note that Ylr^o fr oi"^) ~ 1' ^ should be. The asymptotic limiting value of the correlation 
function obtained above may be understood in the sense of the leading term in a expansion of 
the said quantity (|26p for the fixed case. It should be interesting to obtain successive correction 
terms, and their group theoretic interpretations, to the leading value in the large limit. 



3 The Jacobi matrix of Hahn polynomials 

3.1 General correlation function and special cases 

The method outlined in the beginning of the previous section is clear, and this analysis can in 
principle be made for any set of discrete orthogonal polynomials. So in this section we shall 
consider the more general class of Hahn polynomials Qn{x;a, f3, N) [20,21], characterized by a 
positive integer parameter and two real parameters a and /3 (for orthogonality, one should have 
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a > — 1 and /3 > — 1, or a < —N and /? < —N). The Hahn polynomial of degree n {n = 0,1, . . . , N) 
in the variable x is defined by: 

-n,n + a + f5 + 1, —x 



a + l,-N 



Qn{x) = Qn{x;a,(3,N) = 3F2 
The orthogonality relation reads: 

N 

^w(x)Qn{x)Qra{x) = dn^mn, 

where 



;i 



x=0 



(56) 



(57) 



w{x) 



a + x\fN + /3-x 
N -X 



(x = 0,l,...,iV); 



n\{N -ny.{n + a + (3 + 1)n+i{(3 + l)n 
Ar!2 (2n + a + /3 + l)(a + l)„ 



We have used the common notation for hypergeometric scries and Pochhammer symbols [23,24], 
like {a)n = a{a + 1) • • • (a + n — 1) for n = 1, 2, . . . and (a)o = 1; (a, 6, • • • )„ = (a)„(6)„ • • • , etc. The 
three-term recurrence relation is given by: 



where 



- xQn{x) = An Qn+l{x) - (^n + C'n) Qn{x) + C„ Qn-l{x), 



_ {n + a + p + l)(n + a + l)(iV - n) _ n(n + a + /3 + iV + l)(n + /3) 



(58) 



(2n + a + /? + l)(2n + a + /? + 2) 
Introducing orthonormal Hahn functions 



{2n + a + P){2n + a + 



(59) 



one has the following result [25]: 

Lemma 2 Let Mq be the tridiagonal {N + 1) x [N + l)-matrix (Jacobi matrix) 

\ 



Mq 



ho 


-Jo 







-Jo 


hi 


-Ji ■■■ 







-Ji 


h2 






— JjV-1 

—Jn-1 hN j 



(60) 



where 



Jji — 



{n + l){n + a + l){n + l3+l){n + a + (5+l){n + a + (5 + N + 2){N -n) 
(2n + a + /? + 2)2(2n + a + /? + l)(2n + a + /? + 3) 



N {a- (3) [{a + 13) {N - 2n) - 2n{n + 1)] 



2 ^ 2(2n + a + /3)(2n + a + /3 + 2) 
and let U be the (N + 1) x (A^ + 1) -matrix with elements Ujk = Qj{k). Then 

UU^ = U^U = I and Mq = UDU^ 

where 

D = diag(0,l,2,...,iV). 



(61) 

(62) 
(63) 
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So for a system corresponding to the quantities ([6T]). the transition amphtude is given by 

N 



fr,sit) = E Hk)Qrik)Qsik)z'' (z = e"^*)- (64) 



k=0 



The purpose is now to compute (fMl) . and then to investigate some special cases. Let us denote the 
summation in (|64p by S{r,s): 



N 

S{r, s) = Y, w{k)Qr{k)Qs{k)z\ (65) 

k=0 

In order to perform this summation, one can use the fohowing product formula for Hahn polyno- 
mials: 

^ {-N-(3)k {-k,r-N,s-N,-r-c,-s-c)^ 
(a + l)fc ^0 {l.-N -l3,-c,-N,-N)m 



c - m, 1 + ^^,N + /3 + 1 - m, -m, -r, -s, c + r - iV, c + s - 
-,a + l,c + l,c + l + r — m,c + l + s — m, + 1 — r — m, + 1 — s — m' 



(66) 

where we have used the abbreviation c = N + l + a + p. This expression can be obtained from the 
product formula for (/-Racah polynomials given in [30, eq. (8.3.1)]: in this formula, first take the 
limit a — > 0, and then take the limit q ^ 1. 

Now we multiply the right hand side of (|66p by w{k)z^ , and sum over k from to N. Changing 
the order of summation (over k and m), the inner sum over k can be performed using the binomial 
theorem. This leads to: 



S{r, s) = (-^ra - sFri-l) 

m=0 

{r - N,s - N,-r - N -a- l3 -1,-s - N - a- P -l)r 



m\{-N, -N - [3,-N -a- (5-l)m 



(67) 



where the ^F-j[—l) has the same parameters as in (f66|) . So together with the factor l/\Jdrds 
in (|64p . (|67p gives us a symmetric and compact formula for the computation of fr,s{t) in the Hahn 
case. 

Let us consider the special case when the sender is at one end of the chain, i.e. s = 0. Then ([6^ 
gives 



S{r,0)=^l^Y.(-^ril-z) 



N-m 



{r - N,-r - N - a- P - l)r 



N\ ' ' ' ml(-N-B)m 

(/? + l)Ar f r - N,-r - N - a - f3 - I z 



m 




(1- 


zf 2F, 


(1- 


zY 2F1 



N\ ^ ' ^ ' \ -N - (5 ' z-1 

{P + 1)n,, ,r ^ fr- N,r + a + l 



■,z\. (68) 
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In the last step, Euler's transformation [23, 24] formula was used. Thus, the transition amplitude 
becomes 



/r,0(t) 



{2r + a + P + l){a + l) 



J {13 + l)ria + P + 2)N{r + a + (3 + 
' r — N,r + a + 1 



1/2 



X (p + 1)^(1 - zY 2F1 



-N -13 



; z 



In particular. 



and 



l/jv,o(t)| 



(a + l,/3 + r 



N 



1/2 



{a + 13 + 2)n{N + a + f3 + l)N, 



(1 



sm 



N 



(69) 



(70) 



(71) 



Xa + f3 + 2)NiN + a + 13 + 1)n , 

An interesting special case is that with the parameters a and /3 equal, because then the magnetic 
field strengths hk in ([3]) are all constant (independent of k), see (f6T]) . For /3 = a, ([7T]) becomes 



\fN,om 



(a + r 



N 



{a + 3/2)N-i{a + N/2 + l/2) 



1/2 



N 



(72) 



Clearly, this is maximal for t = it (plus multiples of 2tt). However, \fN,o{t)\ < 1 for the allowed 
values of a. Only for large a, \fN,o{t)\ approaches 1. So "perfect state transfer" between site 
and site does not take place except for a — > +00. This limiting case does not give rise to a 
new example: for a = f3 ^ +00, the Hahn polynomials reduce to Krawtchouk polynomials with 
p = 1/2, and this was the subject of the previous section. 



3.2 A limiting case 

A classical limit of Hahn polynomials are Meixner polynomials. Putting a = b — 1, P = N^-^, and 
letting N go to +00 yields Meixner polynomials Mn{x;b,c) [20]: 

1 — C I — Ti — T 1 \ 

lim Q„(x;6-1,A , N) = Mn{x;b,c) = 2F1 I ' ;1-- , (73) 

c \ b c J 

satisfying the orthogonality relations 

S2^C'Mra{x;b,c)Mn{x-b,c) = " 5mn, (74) 

where 6 > is a positive parameter and < c < 1. For the recurrence relation, see [20]. 

Again one can wonder which spin chain corresponds to this limit. Taking the appropriate limits 
in (fHTI) . one finds for ([3|): 

^ ^k + c{k + b) ^ ^ ^c{k + l){k + b) . t t 

A:=0 fc=0 

It may be immediately observed that the sm(1, 1) algebra acts as the spectrum generating algebra 
of the above Hamiltonian. We define the su{l, 1) generators and the identity operator as bilinear 
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constructs of the fermionic variables by 



A;=0 

00 

k=0 
00 

aiak+i, 

k=0 



(76) 



k=0 



By virtue of ([2|) the above generators satisfy the su{l, 1) commutation relations: 

[/Co,/C±] = ±/C±, [/C+,/C_] = -2/Co, = where X G {/Co, /C±}. (77) 

The large Hamiltonian (I75p now assumes the form 



1 — c 2 1 — c 



(78) 



The infinite-dimensional lowest weight representation of su(l, 1) reads (see, e.g. [31, eq. (2.2)] 
or [32]) 



/Co 



V^(n + 1) (n + 6) 



,n + 1 , 



, = ^/n {n + b- 1) ^, n - 1^ , 



n + 



(79) 



where the lowest weight | has been made explicit in the notation of the state vector. Following [29] 
the BCH-factorization of the time evolution operator (jlip may now be easily obtained: 



U{t) 



exp ( ^/c 



1 — cz 



/C4 



{l-c)V-z \ 
1-cz J 



2/Co 



exp ( ^/c /C- 

l — cz 



(80) 



Employing the decomposition (j80p the correlation function (jl2p for the asymptotic limit of the spin 
chain governed by the Hamiltonian (f78j) is readily obtained. As expected, this precisely agrees with 
the appropriate limit that can easily be computed from ([67j) : 



(I-ZY 



(1 - czf+^'-+'^ 



2FI 



{i-zy 



(81) 



As noted in section [23] the above asymptotic limit of the transition amplitude (|64p may be under- 
stood as its leading term in a expansion scheme. Let us consider an example here, say for 6 = 1 
and c = 1/2. At time t = the state with spin up at position s = is "released" over the infinite 
chain; then at time t = tt it is "spread" as follows: 



/r,o(7r) 



1/^8 



(82) 



(verify that X^J^o /r o(^) ~ decays exponentially over the chain, only to return back to 

its original configuration at time t = 271. 
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4 The Jacobi matrix of dual Hahn and Racah polynomials 

4.1 General computation for dual Hahn polynomials 

Dual Hahn polynomials will play a special role. First of all, the energy spectrum is not linear. 
Secondly, under certain conditions they will allow perfect state transfer. 

Dual Hahn polynomials Rn{X{x);j,6,N) [20,21] are characterized by a positive integer param- 
eter N and two real parameters 7 and 5 (for orthogonality, one should have 7 > — 1 and 5 > —1, 
or 7 < —N and 5 < —N). The dual Hahn polynomial is not a polynomial of degree n in x, but of 
degree n (n = 0, 1, . . . , N) in X[x) = x{x + 7 + 5 + 1): 



Rn{X{x)) = Rn{X{xy,^,S,N) 



— n, —X, a; + 7 + 5 + 1 
7 + 1,-iV 



The orthogonality relation reads: 



N 



^w{x)RniX{x))Rm{X{x)) = dnSmn, 



x=0 



where 



w{x) 



i2x + j + 5 + l)ij + lU-N),m 

(-I)^(X + 7 + 5 + l)iV+l(<5 + l)xX\ 

7 + n\/5 + A/"-n^ 
n 7 V N -n 



(x = 0,l,...,7V); 



The three-term recurrence relation is given by: 

X{x)RniXix)) = AnRn+liX{x)) - {An + C„) Rn{X{x)) + i?„_i(A(x)), 

where 

An = {n + -f + l){n- N), Cn = n{n - 6 - N - I). 
Orthonormal dual Hahn functions are defined by 

^/w{x)Rn{X{x)) 



Rn{X{x)) = 



and then one can deduce: 

Lemma 3 Let Mr be the tridiagonal {N + 1) x [N + 1) -matrix (Jacobi matrix) 

( ho -Jo \ 
—Jo hi — Ji 

-Ji ^2 ■•• 



Mr 



—Jn-1 
—Jn-1 h]sf J 



where 



Jn = V(« + 1) (n + 7 + 1) (TV - n) (5 + iV - n), 
/i„ = (n + 7 + l)(Ar - n) + n{5 + N - n + I). 



(83) 



(84) 



(85) 



(86) 



(87) 



(88) 
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and let U be the {N + 1) x (N + \)-matrix with elements Ujk = Rj{X{k)). Then 

UU^ = U^U = I and Mr = UDU^ (89) 

where 

D = diag(eo, ei, e2, . . . , ej^) with = k{k + 7 + 5 + 1). (90) 

Note that due to the appearance of \{x) in (jSSp . the matrix D has the form (j90p . and thus the 
energy eigenvalues of the single fermion Hamiltonian eigenstates are quadratic in k. 

The rest of the analysis is again concerned with the general correlation function. For a system 
corresponding to the quantities (jSSp . this is now given by 



1 ^ 

UAt) = ^==y^w{k)Rr{\{k))R,{\{k))z^^^+^+'^^^ {z = e-^*). (91) 

In this case, one can use the product formula for g-Hahn polynomials [30, Eq. (8.3.3)] and take the 
limit g — > 1 to find: 

n t\(i.WD f-r,-k,k + -f + 6 + l \ f -s, -k,k + j + 6 + 1 

Rr{X{k))R,{X{k)) = 3F2 7 + 1, -AT ' V ' ' V 7 + 1, -A^ ' 

= ( i)k i^ + ^)k ^ {-k,r-N,s-N,j + S + k + lU 

(7 + l)fc^o il,S + l,-N,-NU 

X F ( -m,-r,-s,-5-m \ 

^ ^\-f+l,N+l-r-m,N + l-s-m' J' ^ ' 

Then one obtains, using (I92p and exchanging the order of summation: 

N 



1 {r-N,s-NU J, ( 

^/drds (0 + l)m V 



{r — N,s — N)m ^ I —m,—r,—s,—S — m 

7 + l,A^+l — r — m, A^+1 — s — m' 



m=0 

^ ( .m iiN-my.f ^ (- jV)fc(7 + S + k+ 1)^(7 + S + 2k + l) ...^ 
^ ^ mlNl ^ {k-m)\{-f + 5 + k + l)N+i ' ^ ' 

Due to the appearance of z^^^^'^^^^^\ the inner sum is no longer of hypergeometric type, and it 
cannot be simplified in general. Let us therefore specialize to the case with s = and r = N 
(sending site at one end and receiving site at the other end of the chain) . Then 

^ k\{j + d + k + 1)n+i 

So far, 7 and 6 are free parameters. Let us now require the following condition: 7 + 5 is an 
odd integer number. Then at time t = ir one has = (_i)'i:(fe+7+5+i) = (-l)'^. The 

summation over k in ([94]) can now be performed, since it corresponds to a nearly-poised 3F2{—1) 
(see [24, (IIL25)]). One obtains: 



Clearly, this expression assumes its maximum value for 7 = 5, and in that case it is equal to 1. In 
other words, for 7 = 5 = p + ^, with p an integer, there is perfect state transfer between the sites 
and N at time t = tt. 
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Similarly, when 7 + 5 is of the form 

7 + ^^=^^, P,gGN ((7/0), (96) 
it also follows that at time t = qir 

^fc(fe+7+5+l) ^ ^-iqn[k(k+l) + ^k] ^ g-i7r(2p+l)fc ^ (-1)^. (97) 

So in that case fN,oiQ'^) assumes the same value as given by the right hand side of ([95|) . Thus for 
7 = (5 = ^^1^, one has perfect state transfer from to at time t = qn. This is the situation 
described in [15]. 



4.2 The case of Racah polynomials 

Dual Hahn polynomials are a limiting case of Racah polynomials. Among the discrete orthogonal 
polynomials, Racah polynomials are the most general, but also the most complicated. Their Jacobi 
matrix is not mirror-periodic, so perfect state transfer is not possible [6]. Racah polynomials 
i?„(A(x); a, /?, 7, (5) are polynomials of degree n in the variable X{x) = x{x + ^ + S + 1), and are 
expressed as a 4^3(1) series, where one of the numerator parameters a + 1, (3 + 6+1 or 7 + I should 
be — A^, with N a positive integer [20] . Without loss of generality, let us assume we are in the first 
case with a + 1 = —N. Then, for the weight function to be positive, assume 

7 + l>0, 5 + l>0, I3>j + N. 

The orthogonality relation and the coefficients of the recurrence relation become fairly complicated, 
see [20]. Moreover, the single fermion eigenvalues are of the same form as (j90p . This means that 
the correlation function is given by (j91ll . where w^k) and dn now stand for the weight function and 
squared norm of the Racah polynomials respectively, and where Rn{X{k)) is a Racah polynomial. 
Due to the appearance of 2:^{^+'y+^+'^) ^ the final summation is again no longer of hypergeometric 
type and cannot be simplified in general. Let us therefore not give the general expression, but only 
some special cases. For s = and r = N, the correlation function becomes 

- VWo k\{j + 6 + N + 2,ij + 6 + l)/2), • ^'"^ 

When ^'i^('':+7+'5+i) = (— l)'^, this sum becomes a nearly-poised ^F2{—1) series which can be summed 
using [24, (III. 25)]. In other words, when 7 + 5 is an odd integer number, one finds 



f I ^ {i + l- [3,5 + 1+ I5)n V(7 + l,^ + l)jv .QQ. 

^"'°^"^ = V (^ + 1). • ^''^ 

Note that in the limit (3 — > +00, in which case the Racah polynomials become dual Hahn polyno- 
mials, pUj) indeed becomes 



5 Discussion and conclusion 

In this paper, we have considered linear spin chains with a nearest-neighbour hopping interaction, 
as models for quantum communication. We have considered the time evolution of single fermion 
states in such a spin chain. In particular, if the system is at time t = in a pure state with all 
spins down except one spin up at site s of the chain, we have studied the behaviour of this system 
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at time t by computing the transition amplitude /r,s(^)- The main contribution of the paper is 
to show that one can deduce closed form expressions for this transition amplitude (or correlation 
function) if the interaction matrix of the system is related to the Jacobi matrix of a set of (discrete) 
orthogonal polynomials. 

We have worked out in detail the cases related to Krawtchouk polynomials (section 2), Hahn 
polynomials (section 3) and dual Hahn polynomials (section 4); for the case of Racah polynomials 
we give only some partial result in section 4.2. Experts in orthogonal polynomials might wonder 
why we did not proceed the opposite way, starting from the most general case (Racah polynomials), 
and then obtaining the other cases as certain limits. This approach would work here only for the 
limit from Racah polynomials to dual Hahn polynomials, because of the appearance of 2:''i^+'y+^+'^) 
in (j98p and (j9ip . For the other correlation functions, there is just in the summation part, and 
these need to be treated separately anyway. The Krawtchouk case could have been presented as a 
limit of the Hahn case; however, we felt it was better to start with a simple example first, which 
has moreover some additional interesting properties (such as the group-theoretical interpretation, 
the special case of perfect state transfer, and an interesting limit of its own). The case of Hahn 
and dual Hahn polynomials had to be considered separately because of the different nature of the 
correlation function (z^ as opposed to 2fc('=+7+'5+i)^^ 

For all examples considered here, we have obtained complete or partial results. In the case of 
Krawtchouk polynomials, we obtained a simple closed form expression for /r,s(t) in general. This 
was also the case for a limit consisting of an infinite chain of spins described by the Jacobi matrix 
of Charlier polynomials. We noticed that this example is related to the unitary representation of 
the Heisenberg algebra. In the case of Hahn polynomials, we do obtain a general expression (|67p 
for fr,s{t), though it is still quite complicated. Some special cases have been discussed, as well as 
the limit related to Meixner polynomials (where a simple general expression is obtained). We have 
related this example with the su{l, 1) symmetry algebra of the corresponding Hamiltonian. The 
asymptotic N ^ co limit of the spin chain may be understood in the framework of a expansion 
where the leading terms of the correlation functions of the related Hamiltonians are obtained. In 
the case of dual Hahn polynomials, the general expression for /r,s(i) contains a summation part 
which is not of hypergeometric type, but we do show that it simplifies in special cases. The same 
remarks hold for the case of Racah polynomials. 

It is worth discussing certain related areas where our analysis may find extensions or appli- 
cations. Propagation of entangled states in anisotropic spin chains has been studied in [16]. 
Anisotropic models are characterized by the property that they allow instantaneous creation of pair- 
wise entanglement from the fully polarized ground state. The anisotropy parameter connects [16] 
the isotropic XY model with the quantum Ising model, which, in the N ^ oo limit, undergoes a 
quantum phase transition at a critical value of the coupling constant. In our case it should be of 
interest to understand the dynamics of propagation of entangled states in a spin chain governed 
by an anisotropic variation of the Hamiltonian ([1]), where the coupling constants follow the poly- 
nomial structures considered here. A class of Hamiltonians that do not preserve the total number 
of excited spins was found [33] to dynamically create multipartite entangled states starting from 
an initially uncorrelated state. In this context one may also introduce the one-axis spin squeezing 
interactions in the Hamiltonian [34] that is expected to protect the entangled states against de- 
coherence. The group theoretic method developed in our work may help in describing analytical 
solutions for evolutions of such entangled states in the quantum register. Work towards this is in 
progress. 

We complete our work mentioning another context in which the formalism developed here 
may be relevant. A recent work [35] considers a Jaynes-Cummings-Hubbard (JCH) system that 
describes coupled cavity structures where confined photons are induced to interact via their coupling 
to embedded two-state systems. In particular, a nonuniform "parabolic" coupling between the 
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cavities is assumed in [35] that is identical to our discussions in section [2] regarding the nearest- 
neighbour coupUng guided by the Krawtchouk polynomials. The large N asymptotic limit of the 
correlation function obtained in (j54p may be useful in understanding the mean field results, and 
consequently, the quantum phase transitions for the JCH Hamiltonians. Moreover, our examples 
of Jacobi matrices corresponding to Hahn polynomials and their asymptotic limits, and the dual 
Hahn polynomials may also be considered as pre-engineered couplings between the cavities for JCH 
systems. Our evaluation of the correlation functions for these cases may have some importance 
in developing the theory of JCH systems with inter-cavity couplings subject to these polynomial 
structures. 
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